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This  report  describes  a  computational  Structure-Battery  Design  Tool  (SBDT)  developed  at  the  Naval  Research  Laboratory  for  analyzing  the 
mechanical  and  electrical  performance  of  multifunctional  structure-battery  materials  configured  in  prismatic  beam  geometries.  SBDT  is  imple¬ 
mented  in  Excel  spreadsheet  form  and  is  capable  of  analyzing  composite  designs  with  several  cross-section  geometries  including:  circular- 
annular,  rectangular-annular,  arbitrary-box,  and  multilayers.  Instructions  for  using  the  SBDT  and  an  overview  of  the  calculations  performed 
therein  are  included  below. 
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Excel  Computational  Design  Tool: 
Multifunctional  Structure-Battery  Materials 


Introduction 

Multifunctional  materials  are  “material  systems”  with  more  than  one  primary  function. 
Multifunctional  materials  are  developed  to  achieve  system-level  performance  enhancements;  local 
subsystem  performance  is  a  secondary  concern.  For  example,  the  best  structure-power  multifunctional 
material  may  not  be  the  strongest  or  stiffest,  or  have  the  largest  energy  storage  capacity.  What  it  should 
achieve  is  better  system-level  performance  (e.g,,  flight  time  for  an  unmanned  air  vehicle)  than  can  be 
obtained  with  any  combination  of  unifunctional  structure  and  power-storage  materials. 

The  Structure-Battery  Design  Tool  (SBDT)  is  implemented  in  the  form  of  Excel  spreadsheets  and 
is  capable  of  analyzing  the  mechanical  and  electrical  energy  storage  properties  for  prismatic  composites 
(Figure  1)  with  circular-annular,  rectangular-annular,  or  arbitrary-box  cross-section  configurations  (Figures 
2).  Inputs  to  the  code  include  information  related  to  geometry,  configuration,  material  properties,  etc.,  and 
outputs  include  mechanical  stiffness  and  strength  under  axial,  bending,  torsion,  shear,  and  buckling 
loads,  electrical  storage  capacity,  weight  and  volume  based  energy  and  power  densities,  etc.  (Figure  3). 
Up  to  five  distinct  materials  can  be  considered  with  the  circular-  and  rectangular-annular  cross-sections 
(Figure  2a  and  2b)  and  seven  distinct  materials  with  the  arbitrary-box  cross-section  (Figure  2c).  The 
analysis  of  mechanical  performance  is  based  on  classical  Mechanics  of  Materials  beam  equations 
suitably  extended  for  composites  through  the  use  of  “modulus-weighted"  cross-section  properties  in  the 
stress  and  deformation  equations.  Nominal  electrical  performance  is  calculated  assuming  constant 
voltage  and  currents  such  that  complete  battery  discharge  occurs  in  one  hour  (i.e.,  1 C  discharge  rate). 


Figure  1:  Schematic  of  a  prismatic  structure-battery  material  with  internal  load  vectors  R  and  M 
(defined  on  page  6,  below)  applied  at  the  ends. 


Manuscript  approved  March  19, 2003. 
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The  type  of  geometric  and  material  input  that  is  required  by  each  template  is  essentially  the 
same.  The  geometric  input  defines  each  material’s  geometry  such  as  width,  height,  placement,  surface 
area,  etc.  The  material  input  is  comprised  of  mechanical  properties  such  as  elastic  stiffness,  strength, 
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density,  and  electrical  properties,  i.e.,  electric  charge  storage  capacities  and  nominal  voltages.  The 
combination  of  the  geometric  and  material  input  is  employed  to  calculate  geometric  and  modulus- 
weighted  parameters  such  as  centroids,  cross-sectional  areas,  areal  moments  of  inertia;  mechanical 
parameters  like  stiffness,  maximum  displacements,  yield  strength,  interfacial  stresses;  and  electrical 
parameters  for  parallel  and  series  connections  like:  voltage,  power  and  energy  densities,  etc.  Each  SBDT 
has  two  calculation  pages.  The  first  page,  denoted  as  “Input-Output”,  performs  calculations  on  one 
particular  design.  The  layout  of  this  page  is  shown  in  Figure  3.  The  second  page,  denoted  “Parametric”, 
has  all  of  inputs  and  outputs  from  “Input-Output”  (page  one)  arranged  horizontally  in  a  single  row.  This 
provides  an  easy  method  for  performing  parametric  studies  on  any  of  the  input  parameters. 


Directions  for  Use 

Colors  are  employed  in  the  structure-battery  design  tool  to  distinguish  input  from  output,  and 
geometric  and  material  data  from  structure-power  material  parameters.  The  color  scheme  is  consistent  in 
all  design  templates. 

The  first  page  in  each  template  is  named  “Input-Output”  and  contains  three  boxes  or  windows. 
The  yellow  box  on  the  left  is  Geometric  and  Material  Data  input  and  calculations,  and  the  tan  box  on  the 
right  is  Structure-Power  properties  for  the  composite.  There  is  a  green  box  below  the  other  two  that  is  not 
shown  in  Figure  3;  it  contains  auxiliary  data  employed  in  calculating  output  in  the  other  two  boxes.  Input 
cells  are  rose-colored  and  output  cells  are  pale  blue.  The  second  page  is  named  “Parametric"  and  uses 
the  same  color  scheme  as  the  first  page. 

The  dimensional  units  used  for  each  input  and  output  cell  are  given.  Note  that  all  Inputs  must 
conform  to  these  units.  Explanatory  comments  are  also  provided  at  various  locations  in  the  template. 

Input-Output:  Geometric  and  Material  Data 

The  circular-annular  design  tool  is  denoted  as  C-SBDT,  the  rectangular-annular  design  tool  is 
denoted  as  R-SBDT,  and  the  arbitrary-box  design  tool  is  denoted  as  AB-SBDT.  The  first  inputs  in  the 
yellow  Geometric  and  Material  Data  box  are:  “Radius”  for  C-SBDT;  “Width”  and  “Height”  for  R-SBDT;  and 
“Width”,  “Height”,  and  the  z'  centroid  locations  for  the  fh  material  regions  (i.e.,  z]  values).  As  the  input  is 
entered,  the  blue  columns  to  the  right  will  be  updated  to  show  information  on  inner  and  outer  diameters, 
wall  thicknesses,  and/or  y'  centroid  locations  for  the  fh  material  region  (i.e.,  y\  values).  Specific  data 
entry  instructions  for  solid  and  hollow  configurations  are  given  next  to  the  cross-section  on  Page  1 . 

The  next  input  consists  of  tensile  and  shear  moduli,  E,  and  G; ,  for  each  component  material. 
The  shear  modulus  value  is  automatically  calculated  using  the  isotropic  material  expression: 
G>  =  E//2(1  +  u)  with  o  (Poisson  ratio)  equal  to  0.33.  The  user  can  replace  the  calculated  shear  modulus 
value  with  actual  data  if  so  desired. 
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The  next  input  is  the  reference  moduli:  ER  and  GR .  Numerical  values  can  be  arbitrarily  assigned. 
A  value  of  1 .0  is  commonly  used  or  the  moduli  values  for  one  of  the  component  materials. 

Yield  strengths  in  tension (a0)t  and  shear  (r0)/  follow.  The  shear  strengths  are  automatically 

calculated  using  the  Tresca  criterion:  (r0).  =  (cx0)./2.  The  user  can  replace  these  calculated  shear 
strength  values  with  actual  data  if  so  desired. 

Component  weight  densities,  p,,  are  next.  The  last  box  in  this  column  is  output:  the  overall 
weight  density  of  the  composite. 

The  next  input  is  surface  area  per  unit  beam  length  (i.e.,  lineal  surface  area),  S, .  Each  S,  entry 
is  currently  assigned  the  value  of  the  average  circumference  or  perimeter  of  the  material  region  (e.g., 
*(n+rM)  for  circular  cross-sections).  The  user  can  replace  these  calculated  values  with  actual  data  if 

so  desired.  The  last  two  cells  in  the  column  are  the  total  sum  of  lineal  surface  areas,  S  ,  and  the  lineal 
surface  area  of  the  outer  layer,  S0D.  Surface  area  quantities  are  used  for  electrical  performance 
calculations  on  material  regions  comprised  of  multiple  layers  (e.g.,  multiple  battery  material  layers). 

The  next  two  blocks  are  output.  The  first  consists  of  the  cross-sectional  area  of  each  material 
component,  A( ,  and  total  material  cross-sectional  area  of  the  composite,  A .  The  second  block  gives  the 

material  component  weight  per  unit  beam  length  (i.e.,  lineal  weight  densities),  wLi,  and  the  total  beam 
weight  per  unit  beam  length,  wL . 

Modulus-weighted  (MW)  section  properties  are  denoted  by  the  “  *  ”  superscript.  This  includes 
parameters  like:  material  area,  A' ;  cross-section  location  of  the  modulus  weighted  centroid,  y'*  and  z'*; 

2nd  moments  of  area,  / ^ , /  *z ,  and  ;  the  two  principal  2nd  moments  of  area,  l{2 ,  and  the  corresponding 
angle  between  the  1-2  principal  axes  relative  to  the  y-z  MW  centroidal  axes,  6P ;  polar  moment  of 
inertia,  K‘\ and  1st  areal  moments  of  inertia,  Q'  and  Q'  about  the  MW  cross-section  axes. 

The  last  two  input  blocks  are  electrical  data:  electrical  charge  storage  capacity  per  unit  area,  C, , 
and  nominal  voltage,  V; ,  for  each  material  component.  The  remaining  blocks  are  electrical  property 
outputs  calculated  for  the  fh  material  regions  including:  electrical  charge  storage  capacity  per  unit  length 
(lineal  capacity),  cLI ;  power  capacity  per  unit  length  (lineal  power  density),  pLi ;  power  per  unit  volume 

(power  density),  P, ;  power  per  unit  weight  (specific  power),  p, ;  electrical  energy  per  unit  length  (lineal 
energy  density),  eL/ ;  electrical  energy  per  unit  volume  (energy  density),  E, ;  and  electrical  energy  per 
unit  weight  (specific  energy),  e, .  All  of  the  electrical  property  calculations  assume  a  1 C  current  draw  rate 

(i.e.,  current,  in  amperes,  numerically  equal  to  the  total  charge  capacity  in  ampere-hours  [Linden  and 
Reddy,  2001]). 
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Structure-Power  Material  Parameters 

The  first  block  in  the  tan  colored  box  contains  elastic  stiffness  estimates  for  the  composite  beam. 
The  entries  in  C-SBDT  include  axial,  bending,  and  torsional  stiffnesses.  Bending  stiffness  values  about 
the  y  and  z  axes  are  included  with  R-SBDT.  Bending  stiffnesses  about  y,  z,  and  the  two  principal 

axes  are  included  with  AB-SBDT.  Note  that  the  modulus-weighted  y-z  axes  are  principal  for  the  circular 
and  rectangular  cross-sections. 

The  block  below  gives  elastic  deformations  for  beams  of  unit  length  with  unit  applied  loads.  The 
entries  in  C-SBDT  include  axial  and  bending  displacements,  and  bending  and  torsional  rotations.  R-SBDT 
includes  bending  displacements  and  rotations  about  the  y  and  z  axes,  and  AB-SBDT  includes  bending 
displacements  and  rotations  about  the  y ,  z ,  and  two  principal  axes. 

The  block  to  the  right  gives  the  electrical  performance  characteristics  for  the  composite  beam  in 
the  parallel  and  serial  configurations.  This  includes  quantities  like:  voltage,  lineal  and  areal  charge 
storage  capacities,  and  lineal,  volumetric,  and  weight-based  energy  and  power  densities. 

Yield  data  are  given  in  the  next  row  of  blocks.  This  includes  yield  values  for  axial  force,  bending 
moment,  torque,  and  transverse  shear  for  each  material  and  for  the  composite  as  a  whole.  R-SBDT 
includes  two  bending  moment  and  transverse  shear  yield  loads  corresponding  to  y  and  z  loading.  AB- 

SBDT  includes  four  bending  moment  yield  loads  corresponding  to  y ,  z ,  and  principal  axis  bending.  The 

last  entry  in  each  block  is  the  composite  beam  yield  load,  and  it  is  taken  as  the  minimum  of  all  component 
material  yield  loads. 

The  next  “block”  gives  the  Euler  buckling  forces  for  composite  beams  of  unit  length.  Values  are 
computed  for  buckling  about  the  principal  2nd  moment  of  area  axes. 

The  last  row  of  blocks  displays  the  normal  and  shear  stress  discontinuities  at  the  interface 
between  component  materials  with  unit  internal  load  values  (i.e,,  axial  force,  bending  moment(s)  and 
torque).  The  last  entry  is  the  maximum  stress  discontinuity  value  for  that  particular  loading. 


Parametric  Page 

The  entries  on  the  Parametric  page  are  exactly  the  same  as  those  on  the  Input-Output  page 
except  for  their  physical  layout.  The  cells  on  the  Parametric  page  are  arranged  in  a  single  horizontal  row 
of  the  spreadsheet.  Parametric  studies  can  be  performed  by  copying  a  “fiiled-in"  row  (i.e.,  inputs  and 
subsequent  outputs)  to  one  or  more  empty  rows  followed  by  modifications  to  one  or  more  of  the  column 
parameters.  After  performing  the  copy-paste  and  parameter  modification  steps,  copy  the  entire  page  and 
then  use  the  “Paste-Values  Only”  command  to  copy  the  only  the  numerical  values  onto  another  page  for 
use  in  plotting.  The  "Paste-Values  Only”  command  allows  one  to  delete  useless  columns  without  affecting 
the  other  numbers.  On  the  other  hand,  if  changes  are  made  on  the  Parametric  page,  then  the  page-copy- 
for-plotting  steps  will  have  to  be  re-performed. 
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Mechanical  Analysis  of  Prismatic  Structure-Battery  Materials 

The  use  of  modulus-weighted  cross-section  properties  provides  a  simple  method  for  taking  into 
account  the  composite  nature  of  prismatic  structure-battery  materials  in  the  analysis  of  mechanical 
performance  [Allen  and  Haisler,  1985].  Standard  “Mechanics  of  Materials”  beam  equation  forms  are 
retained  with  certain  cross-sectional  parameters  replaced  by  modulus-weighted  quantities. 

The  materials  and  cross-section  are  assumed  to  be  constant  along  the  length  of  the  member 
(prismatic  beam  assumption),  taken  as  the  x-axis.  Six  types  of  internal  loads  are  possible.  They  are 
assumed  to  be  constant  along  the  length  (i.e.,  no  x  dependence)  and  to  pass  or  act  through  the  cross- 

section’s  modulus-weighted  centroid.  They  are  represented,  in  Figure  1,  by  the  vectors  R  and  M : 

R  =  Px  +  Vyy  +  Vzz  (1) 

— > 

M  =  7x  +  Mxyy  +  Mxzz  (2) 

P  is  the  axial  force;  Vy  and  Vz  are  the  y  and  z  transverse  shear  forces,  respectively;  T  is  the  torque; 
and  Mxy  and  Mxz  are  the  y  and  z  bending  moments,  respectively.  The  internal  load  sign  convention  is: 

Positive  internal  loads:  positive  face  and  positive  coordinate  directions  or  negative  face  and 
negative  coordinate  directions. 

Negative  internal  loads:  positive  face  and  negative  coordinate  directions  or  negative  face  and 
positive  coordinate  directions. 


Modulus-Weighted  Cross-Section  Property  Calculations 

The  modulus-weighted  cross-section  properties,  for  discrete  material  distributions,  are  calculated 
as  shown  below  [Allen  and  Haisler,  1985].  The  y'-z' coordinate  axes  can  be  arbitrarily  located  in  the 

plane  of  the  cross-section.  The  primary  coordinate  system  is  the  y-z  modulus-weighted  centroidal  axes. 

Modulus-Weighted  Area: 


/=1  CR 


Modulus-Weighted  Centroid  Location: 


— 1  A  E  1  A  £: 

y  :=  and  2 

A  /= ItR  A  m  tR 


Modulus-Weighted  2nd  Area  Moments  of  Inertia: 


C  +Z?Ai]  '  <n  +yfAl]  -  C  ■=t§L{(l,o)i  +y<Z,Ai] 

/= 1  CR  /=1  CR  /=1  CR 

f 


1 1.2  := 


/„+/ 


yy  ’  '  zz 


/»-/ 


yy  'zz 


+  (lyz)  and  26>p  =  tan  1 
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yz 


\lyy  I  zz  ) 


(3) 


(4) 


(5) 


(6) 
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(7) 


Modulus-Weighted  Polar  Moment  of  Inertia  (parallel  torsional  springs); 

SfKKo), 

In  Equations  (3)-(7),  n  is  the  number  of  materials  in  the  cross-section.  E;  and  G,  are  the  tensile 
and  shear  modulus  for  the  fh  material.  A ,  is  the  cross-sectional  area  of  the  material,  y]  and  zj  are  the 
coordinates  of  the  P1  material  centroid  relative  to  the  y'-z'  axes,  (/s0).,  and  (lyz0)i  are  the 

second  moments  of  area  of  the  ih  material  about  the  fh  material  centroidal  axes,  y,  and  z,  are  the 
coordinates  of  the  /  material  centroid  relative  to  the  modulus-weighted  y-z  centroidal  axes.  /12  are  the 
principal  2  d  moments  of  area  and  6p  is  the  principal  orientation  angle  relative  to  the  y-z  centroidal  axes. 
(K0),  is  the  “effective"  polar  moment  of  inertia  of  the  fh  material  area  relative  to  the  fh  material  centroidal 
axes.  Expressions  for  (K0).  pertaining  to  rectangular  cross-sections  are  given  below.  Finally,  Ep  and 

Gr  are  the  reference  tensile  and  shear  moduli;  their  value  is  arbitrary  and  any  convenient  value  can  be 
selected  for  use. 

Note  that  1^  =  0  for  the  circular-  and  rectangular-annular  cross-sections,  for  doubly  symmetric 

cross-sections,  and  for  cross-sections  with  the  modulus-weighted  coordinate  axes  in  the  principal 

orientation.  Great  simplifications  occur  in  the  stress  and  deformation  expressions  when  this  cross-product 
term  is  zero. 


Stress  Calculations 

The  normal  stress  in  the  fh  material  due  to  internal  axial  force  and/or  bending  moment  is  given  by: 

P  MJl+Mjl  MJl+MJ 


xz 1  yz 


(CC-C1) 


(8) 


where  “y”  and  "z”  are  positive  or  negative  distances  from  the  modulus  weighted  centroid  fo  the  point  of 
interest. 


The  shear  stress  in  the  fh  material  due  to  an  internal  torque  is  given  by: 

G,  T xr 

r,=G^r 

for  circular-annular  cross-sections  where  V  is  the  distance  from  the  modulus  weighted  centroid  to  the 
radial  point  of  interest.  The  maximum  torsional  shear  stress  for  solid  rectangular  sections  is  given  by: 


_G**,K*»,T(3b  +  1.8h) 
Gr  K'  b2h2 


(10) 


It  occurs  on  the  perimeter  at  the  midpoint  of  the  longest  side.  The  geometry  is  shown  in  Figure  4a  with  the 
“ b *  dimension  denoting  the  longer  side. 
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The  torsional  shear  stress  for  hollow  rectangular  sections  is  given  by: 

(r.)  =5^ _ I _ 

'  Gr  K-  2tu(b-t,)(h-t,) 

and  is  assumed  constant  through  the  wall  (thin-walled  assumption).  The  geometry  is  shown  in  Figure  4b. 
Expressions  for  the  “effective”  polar  moment  of  inertias  are  given  by: 


K  solid  ~  bh 


0.21"  1- 


nMt>-tb)2(h-th)2 

btb+hth-tb2-th2 


(13) 


Eqs.  (10)  and  (12)  are  obtained  from  [Roark  and  Young,  1975]  and  Eqs.  (11)  and  (13)  from  torsion 
analysis  of  arbitrary  shaped  thin-walled  tubes  [Allen  and  Haisler,  1 985], 


(a)  (b) 

Figure  4:  Solid  (a)  and  hollow  (b)  rectangular  cross-section  geometries  and  variables. 


Shear  stresses,  rxy  and  ,  due  to  transverse  shear  forces  are  given  by: 

.  _kc^,c)q;  dnd  .  K(;*v;)s; 

(cc-/;2)  *.  t,  ( 

where  Q’y  and  Q‘  are  the  modulus-weighted  first  moments  of  area  about  y-z  axes: 

a"*  (ie 

f= 1  /=1  *—  R 

with  (Qy).  :=  Z/Aj  and  (Qz).  :=  y,A, ,  the  first  moments  of  area  of  the  fh  material  about  the  modulus- 


weighted  centroidal  axes  (i.e.,  y  and  z).  Also,  ty  and  tz  are  the  material  thickness  (total)  along  y  or  z 
cross-section  “cuts”  where  the  shear  stresses  are  computed. 
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Yield  Failure  Loads 

Failure  by  yielding  is  assumed  to  occur  when  any  one  of  the  component  materials  yields.  Internal 
load  values,  considered  one  at  a  time ,  for  which  yielding  occurs  are  calculated  using  the  stress 

equations,  Eqs.  (8)-(10),  (1 1),  and  (14),  and  tensile  and  shear  yield  strength  input  values.  The  equations 
are  detailed  below: 


C-SBDT : 


R-SBDT: 


(po ),  =-jjr-A'(cr0)i  and  P0  :=min(P0). 

e  r 

iMo)i=~(  -  T  (aQ).  and  M0:=mm(M0)l 

cf  V  max  }  i  1 

Q 

^0^  =  G5'7r~r(r°)''  and  To-=mm(T0)l 

V  max  //  1 

(v° )'  =  t~Q^(T°)i  and  Vo^minfVo), 

OM,  =|5-77JVK)/  and  Mxy0:=min(Mxy0) 

e  r 

(M*zo),  (°~o)f  and  Mxz0  :=min(W„0). 

max  /  /  1 


AB-SBDT: 


(7o)so/W  G^KS0(a.(3b  +  1.8/j)^°^«  '  (To)i-^Jj^2tP*(b-h)(h-th)(r 0), 

T0:=min{(T0)soM,(T0)(} 
t  I* 

(Vy°)i  =^f-(r°)/  and  ^yo  :=min(\/r0), 

(^o),.=^l(r0)f  and  t/z0  :=  min(\/z0). 

r*  OM,  =§“77z  l"  ”,  *  \ 7 \(ffo)/  and  Mxyo:=min(M,y0) 

f  ^^max/^zz  \y  max ) \* yzj 

(^xzo)i  -  E  77-  t  .  \  .  Jo'o),  and  M„0  :=min(M„0) 

1  \\4m ax//fyz  ^/niax//'yyj 

(Mio),=|g-p^-j  (o-0),  and  M10:=min(M10). 

e  r 

(M™)i=lT7^r  (°o),  and  M20  :=min(/W20) 
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(T0).=— ■ — — — — -(r0).  and  T0  :=  min{(T0).} 
v  0>l  G,  K,  (3ft +  1. 8ft) v  h  '  U  o;,/ 

(U  >(^o),  and 


yy 


(^o)j=  ^  (VfiT_^J(r°)(  and  \Zz0:=min(\/z0). 

Wy  /  ZZ 


(30) 

(31) 

(32) 


Buckling  Failure  Forces 

Critical  buckling  forces  for  the  principal  planes  (xy  and  xz  planes  for  the  circular  and  rectangular 
cross-sections)  are  calculated  using  the  Euler  buckling  equation  with  pinned  end  conditions  and  unit 
beam  length  (i.e.,  L  =  1  cm ): 


(Pcr),= 


*2eri; 

L2 


and  (PCr)2  = 


7t 


lERi2 


(33) 


The  minimum  buckling  load,  for  a  given  cross-section  configuration,  corresponds  with  the  minimum 
principal  2nd  areal  moment  value. 


Interphase  Stress  Discontinuities 

The  difference  in  the  normal  and  the  torsional  shear  stresses  at  the  interface  between  the 
component  materials  is  calculated  for  various  unit  internal  loading  cases. 


Deformation  Calculations 


Displacements  and  rotations  for  the  composite  beams  are  calculated  assuming  that  the  internal 
loads  are  constant  along  the  length  (i.e.,  no  x  dependence): 


do, , 

dx 


ErA 


=>  u„  =■ 


PL 

ErA' 


d2v„ 


M  xz^  yy  +  ^  Xyl  yz 


dx2 

M/x-/;2) 

d2w0  _ 

M  xyf  zz  +  M  XZI  yz 

dx2 

OS,  _  T 

dx  GrK  ' 

ddy  _ 

Myylzz  +MXZlyz 

W„ 


=>  ex  = 


( Mjyy+M^ ;)L2 

=  2er(/;/;-/;2) 

'  2er(/;/;z-/;2) 

TL 


GrK 


dx  er(/;c-/;2) 


(Mxylzz  +  )/- 

e„(/;/;-/;2) 


(34) 

(35) 

(36) 

(37) 

(38) 
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(39) 


+  Mxyl^yz  ,  xzC  ^  ^  d  yz)L 

dx  ~er(cc-c2)  2"  er(/;/;-/;2) 

In  the  above  equations,  u0)v0,  and  w0are  the  x,  y,  z  displacements,  and  0x,0y,  and  02  are  the  x,  y  and 
z-axis  rotations  of  the  modulus-weighted  centroid. 


Elastic  Stiffness  Calculations 

Elastic  stiffness  values  for  the  various  internal  loads  are  calculated  using  Eqs.  (34)  -  (39): 


Axial: 

p  =  kau0 

=>  := 

era* 

L 

(40) 

Torsion: 

T  ~kt0x 

=>  kt  ■= 

GrK* 

L 

(41) 

Bending-Symmetric: 

M  gy 

=  (k»)yy0y  =* 

(k»)yy  := 

(4C- 

Cl 

■V) 

(42) 

Mxz 

=  {kb)Z20z  => 

<M.  := 

Er  (l yylZz  ~ 

CL 

V) 

(43) 

Antisymmetric:  M„  = 

(*>)**' 

and  Mxx=(kb)zy 

^  (^*) 

yZ=(kt)zy 

.  E,(I'X<2) 

Cl 

(44) 

The  values  listed  in  the  Elastic  Stiffnesses  block  are  computed  assuming  unit  length  (i.e,,  L  =  1 
cm).  Also,  note  that  (kb)yy=(kb)a  and  (k,,)^  =  (kb)zy  =*>  for  circular  configurations  and  that 

(**)„=(**)„=«>  for  rectangular  configurations. 


Elastic  Deflection  Calculations 

Elastic  deformations  are  calculated  for  various  internal  loadings  assuming  unit  load  and  unit 
length  values  (e.g.,  P  =  1  N  and  L  =  1  cm ): 

U°  =  TT  =*  (t/0)0  :=  =  _  -r  (45 


MJ-  .  MxyL 


'”'2(*:l+2(4  "  and  w-=ii4^7) 


-MxyL  -MXZL 
*(*.). +2(*.)„ 


2Ei(rJi-V)  and 
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*y 


dv  =; 


0  =  + 


-Mxz 

^  (^y)s 

-c 

and  (0y)g 

(49) 

(Myz 

er(/w/;z-C2) 

Mxy 

=*  (**). 

/; 

and  (0z)a 

(50) 

(*.)„ 

MC'.-C’) 

Summarizing: 

<?a  is  the  x  deflection  for  P  =  1  /V  and  L  =  1  cm 

8y  or  (<Sy)s  is  the  “symmetric”  part  of  the  y  deflection  for  Mxz  =  1  N-cm  and  L  =  1  cm 
(<Sy)fi  is  the  “antisymmetric”  part  of  the  y  deflection  for  /Wxy  =  1  /V-cm  and  L  =  1  cm 
S2  or  (<5z)s  is  the  “symmetric”  part  of  the  z  deflection  for  Mxy  =  1  N-cm  and  L  =  1  cm 
(Sz)g  is  the  “antisymmetric”  part  of  the  z  deflection  for  M„  =  1  N-cm  and  L  =  1  cm 
6X  is  the  x-axis  rotation  for  T  =  1  AZ-c/t?  and  L  =  1  cm 

6y  or  (/9y)s  is  the  “symmetric”  part  of  the  y-axis  rotation  for  Mxy  =  1  N-cm  and  L  =  1  cm 
(dy)e  is  the  “antisymmetric”  part  of  the  y-axis  rotation  for  Mxz  =  1  N-cm  and  L  =  1  cm 
or  (0z)s  is  the  “symmetric"  part  of  the  z-axis  rotation  for  Mxz  =  1  N-cm  and  L  =  1  cm 
(6>z)a  is  the  “antisymmetric”  part  of  the  z-axis  rotation  for  Mxy  =  1  N-cm  and  L  =  1  cm 
Cross-deflections  and  rotations  (i.e. ,  symmetric  and  antisymmetric  parts)  occur  when  *  0 . 


Weight  Density  Calculations 

The  weight  density  of  the  composite  beam  is  defined  as: 

5>A 

P=J~A  <51> 

where  p{  is  the  weight  density  of  the  fh  material,  At  is  the  cross-section  area  of  the  iih  material  region,  and 
A  is  total  cross-section  area.  The  weight  density  per  unit  length  (lineal  density),  is  defined  as: 

wt:=PA  (52) 


Battery  Performance  Calculations 

The  electrical  charge  storage  capacity  per  unit  area  (areal  charge  capacity),  C, ,  and  the  nominal 
voltage,  V , ,  of  the  fh  material  are  supplied  as  input  data  in  the  yellow  Geometric  and  Material  Data  box. 
Areal  charge  storage  capacity  is  a  common  parameter  used  in  rating  the  performance  of  battery  cell 
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materials  supplied  in  thin-sheet  form.  The  “theoretical”  charge  storage  capacity,  or  the  total  amount  of 
electrical  charge  transferred  in  complete  reaction  of  active  battery  material,  can  be  calculated  knowing  the 
amount  of  active  material  present  in  the  cell  and  the  chemical  and  electron  transfer  reactions  involved. 
Charge  storage  capacity  per  unit  length  (lineal  charge  capacity)  is  output  in  an  adjacent  block: 

Lineal  Charge  Capacity;  Cu  :=  C,S,  (5; 


S,  is  the  surface  area  of  the  fh  material  per  unit  beam  length  (lineal  surface  area). 

The  next  row  of  blocks  gives  the  power  capacity  per  unit  length  (lineal  power  density),  the  power 
per  unit  volume  (power  density),  and  the  power  per  unit  weight  (specific  power)  for  each  materia! 
component.  The  calculated  power  quantities  assume  that  discharge  occurs  at  the  nominal  voltage  and  at 
the  1C  current  draw  rate  (i.e,,  current  in  A  numerically  equal  to  charge  storage  capacity  in  A-h)  (Linden  & 
Reddy,  2001).  At  the  1C  discharge  rate,  the  battery  energy  is  completely  drained  in  1  hour. 

Lineal  Power  Density:  D  •-  c^’ 

u  '  1  hr  < 54> 

Power  Density:  p  ._£u. 

'  -  A,  <55) 


Specific  Power:  p,:=  —  ^ 

The  last  row  of  blocks  gives  the  stored  energy  per  unit  length  (lineal  energy  density),  the  stored 
energy  per  unit  volume  (energy  density),  and  the  stored  energy  per  unit  weight  (specific  energy)  for  each 
material  component.  Again,  the  calculated  quantities  assume  that  discharge  occurs  at  the  nominal 
voltage  and  at  the  1  C  current  draw  rate: 


Lineal  Energy  Density:  e0  :=  pu  x  1  hr 

Energy  Density:  E,  :=  P,  xl  hr 

Specific  Energy:  ef:=p,.x1/jr 

Parameters  in  the  Electrical  Performance  Characteristics  block  located 


(57) 

(58) 

(59) 

in  the  tan-colored 


Composite  Structure-Battery  Properties  box  are  calculated  for  both  parallel  and  serial  configurations 
and  include:  electrical  charge  storage  capacity  per  unit  length  (lineal  charge  capacity),  voltage,  1  C 
current  capacity  per  unit  length  (lineal  current  density),  power  per  unit  length  (lineal  power  density),  power 
per  unit  enclosed  (C-  and  R-SBDT)  volume  (power  density),  power  per  unit  weight  (specific  power), 
stored  energy  per  unit  length  (lineal  energy  density),  energy  per  unit  enclosed  (C-  and  R-SBDT)  volume 
(energy  density),  and  energy  per  unit  weight  (specific  energy)  for  the  composite  beam.  The  equations 
used  in  the  calculations  are  detailed  below: 
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Parallel 

Series 

Lineal  Charge  Capacity: 

CL  =2>u 

i 

CL  =  CLi  =CL2  =  CL  3  =  ••• 

(60) 

Voltage: 

11 

11 

> 

11 

II 

II 

(61) 

Lineal  1C  Current  Density: 

i  =  ^ 

L  1  hr 

,l=a_ 

L  1  hr 

(62) 

Lineal  Power  Density: 

cLV 

c.V 

Pl  =  1  hr 

(63) 

Power  Density: 

p  =  ^ 

A0 

(64) 

Specific  Power: 

-O 

II 

It3 

1  r- 

P 

(65) 

wL 

wL 

Lineal  Energy  Density: 

eL  =  CLV 

et=cLV 

(66) 

Energy  Density: 

E  =  — 

E  =  — 

(67) 

Parallel 

Series 

Specific  Energy: 

eL 

e  =  — 

eL 

e  =  — 

(68) 

wL  wL 

The  cL  (Series)  quantity  in  Eq.  (60)  is  assigned  the  first  nonzero  cu  value  and  the  V  (Parallel) 


quantity  in  Eq.  (61 )  is  assigned  the  first  nonzero  V,  value.  The  quantity,  A0 ,  denotes  the  cross-section 
area  defined  by  the  outermost  perimeter,  for  the  circular  and  rectangular  cross-sections,  and  total 
material  area  for  arbitrary-box  sections. 
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Nomenclature  Table 


A 

A 

A 

A * 

b 


CLI 


CL 

C, 


e 

«i 

eL 


e 

E, 


Ei .  G/ 

Er.Gr 

G  solid 


h 

k 


yyyo)i ■  (lno)j  >  (lyz o), 


Total  material-filled  cross-section  area  of  the  composite  [cm2] 

Cross-sectional  area  of  fh  material  region  [cm2] 

Projected  cross-section  area  of  the  outer  perimeter  [an2] 

Modulus-weighted  cross-sectional  area  [cm2] 

Width  of  a  rectangular  solid  component  [cm] 

Charge  storage  capacity  per  unit  length  (lineal  capacity)  of  /*  material  region  [Ah/cm] 
Charge  storage  capacity  per  unit  length  (lineal  charge  capacity)  of  composite  [Ah/cm] 
Charge  storage  capacity  per  unit  area  (areal  charge  capacity)  of/*  material  [Ah/cm2] 
Stored  energy  per  unit  weight  (specific  energy)  of  composite  [Wh/Kg] 

Stored  electrical  energy  per  unit  weight  (specific  energy)  of/*  material  region  [Wh/Kg] 
Stored  energy  per  unit  length  (lineal  energy  density)  of  composite  [Wh/cm] 

Stored  electrical  energy  per  unit  length  (lineal  energy  density)  of  /*  material  [Wh/cm] 
Stored  energy  per  unit  volume  (energy  density)  of  composite  [Wh/cm3] 

Stored  electrical  energy  per  unit  volume  (energy  density)  of  /*  material  region  [Wh/cm3] 
Tensile  and  shear  modulus  of  (h  material  region  [MPa] 

Reference  tensile  and  shear  modulus  [MPa] 

Shear  modulus  of  rectangular  solid  component  [MPa] 

Height  of  a  rectangular  solid  component  [cm] 

1  C  current  capacity  of  the  composite  per  unit  length  (lineal  current  capacity)  [A/cm] 

2  areal  moments  of  the  /*  material  region  about  their  centroidal  axes  [cm4] 


r  r 
'1  ■  '2 

A 

(**)„■(**)« 


Principal  modulus-weighted  2nd  moments  of  area  [cm4] 
Modulus-weighted  2nd  moments  of  area  [cm4] 

Axial  stiffness  of  composite  [N/cm] 

Torsional  stiffness  of  composite  [N-cm] 

Symmetric  bending  stiffness  components  [N-cm] 


Antisymmetric  bending  stiffness  components  [N-cm] 


(*.), 
k  solid 

K" 

M  Xy  ,  M  ^ 

M xzO  i  M xyQ 
(^xzo)/!  (^xyo)j 


p 

Pi 

Pi 


“Effective"  polar  moment  of  inertia  of  /*  material  region  about  their  centroidal  axes  [cm4] 
“Effective”  polar  moment  of  a  rectangular  solid  component  [cm4] 

Modulus-weighted  “effective”  polar  moment  of  inertia  [cm4] 

Internal  bending  moments  about  y  and  z  axes  [N-cm] 

Bending  moment  components  that  cause  yielding  in  composite  [N-cm] 

Bending  moment  components  that  cause  yielding  in  /*  material  region  [N-cm] 

Power  per  unit  weight  (specific  power)  of  composite  [W/kg] 

Power  per  unit  weight  (specific  power)  of  /*  material  region  [W/kg] 

Power  per  unit  length  (lineal  power  density)  of  composite  [W/cm] 
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Pu 

P 

P 

P, 

P0 

(Po), 

(P»)v  ( Pcr)2 
(Qr),-(Q.), 

q;.q; 

r 

S 

$, 

Sqd 

tb 

th 

ty •  fz 
T 

T0 

(To), 

U0,  V0,  w0 
V 
V, 

vy,  V, 

Vy0,V20 

M,  ■  (v»), 

WL 
W  Li 

x,  y,  z 
x',  y\  z' 

y.  * 

7/  -  ?/ 
y;> 


Power  capacity  per  unit  length  (lineal  power  density)  of  fh  material  region  [W/cm] 

Power  per  unit  volume  (power  density)  of  composite  [W/cm3] 

Internal  axial  force  [N] 

Power  per  unit  volume  (power  density)  of  !h  material  region  [W/cm3] 

Axial  force  that  causes  yield  in  composite  [N] 

Axial  force  that  causes  yield  in  fh  material  region  [N] 

Euler  buckling  forces  for  bending  in  the  principal  planes  [N] 

1st  areal  moment  of  fh  material  region  [cm3] 

Modulus-weighted  1st  areal  moments  [cm3] 

Radial  coordinate  from  modulus  weighted  centroid  to  point  of  interest  [cm] 

Total  sum  of  lineal  surface  areas  [cm2/cm] 

Surface  area  per  unit  beam  length  (lineal  surface  area)  of  fh  material  region  [cm2/cm] 
Surface  area  of  the  outer  layer  [cm2/cm] 

Thickness  of  rectangular  cross-section  of  fh  material  region  in  width  direction  [cm] 
Thickness  of  rectangular  cross-section  of  fh  material  region  in  height  direction  [cm] 
Material  thickness  in  the  y  and  z  directions  [cm] 

Internal  torque  [N-cm] 

Torque  that  causes  yield  in  composite  [N-cm] 

Torque  that  causes  yield  in  fh  material  region  [N-cm] 

Displacement  components  of  the  MW  centroid  [cm] 

Nominal  voltage  of  the  composite  [Volts] 

Nominal  voltage  of  fh  material  region  [Volts] 

Internal  transverse  shear  forces  in  y  and  z  directions  [N] 

Transverse  shear  force  components  that  cause  yielding  in  composite  [N] 

Transverse  shear  force  components  that  cause  yielding  in  fh  material  region  [N] 

Weight  per  unit  length  (lineal  weight  density)  of  composite  [g/cm] 

Weight  per  unit  beam  length  (lineal  weight  densities)  of  fh  material  region  [g/cm] 

Unit  vectors  in  the  modulus  weighted  centroidal  coordinate  system 

Units  vectors  in  an  arbitrary  coordinate  system 

Coordinates  from  modulus  weighted  centroid  to  point  of  interest  [cm] 

Centroidal  coordinates  of  the  fh  material  region  relative  to  the  modulus-weighted  axes  [cm] 
Centroidal  coordinates  of  the  fh  material  region  relative  to  the  y'-z'  axes  [cm] 


$a »  8y  .  Sz 
u 
P 
Pi 


Displacement  components  of  a  point  located  at  the  MW  centroid  [cm] 
Poisson’s  ratio  [1] 

Weight  density  of  composite  [g/cm3] 

Weight  density  of  fh  material  region  [g/cm3] 

Tensile  yield  strength  of  fh  materia!  region  [Pa] 
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*■/ 

^"max 
^xy  * 

(ro)j 

(r/)w 

*x-  '  @z 

W. 

(*)a 


Torsional  stress  in  circular  annular  cross-section  of  material  region  [Pa] 
Maximum  torsional  stress  in  rectangular  solid  component  [Pa] 

Transverse  shear  stress  components  [Pa] 

Shear  yield  strength  of  f1  material  region  [Pa] 

Torsional  stress  in  rectangular  cross-section  of  f  material  region  [Pa] 
Orientation  of  the  principal  axes  relative  to  the  y-z  axes  [deg] 

Rotation  components  of  a  point  located  at  the  MW  centroid  [deg] 
Symmetric  component  of  the  quantity  in  parentheses 
Antisymmetric  component  of  the  quantity  in  parentheses 
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